Back

RNA Biology

Informa UK Limited

Preprints posted in the last 90 days, ranked by how well they match RNA Biology's content profile, based on 78 papers previously published here. The average preprint has a 0.06% match score for this journal, so anything above that is already an above-average fit.

1
Identification of rare pseudouridylated microRNA by comprehensive small RNA bisulfite sequencing of mouse and human tissues

Fagre, C.; Gilbert, W. V.

2026-05-15 molecular biology 10.64898/2026.05.14.725264 medRxiv
Top 0.1%
23.1%
Show abstract

Pseudouridine ({Psi}) is an important post-transcriptional modification of many noncoding RNAs that is under-characterized in microRNA (miRNA) due to historical limitations in pseudouridine mapping methods. {Psi} modification stabilizes RNA duplex structures and could therefore play an important role in miRNA target binding and repression. To investigate the extent to which mammalian miRNAs are modified with {Psi}, we profiled the modification landscape of short (<30 nt) RNA in human cells and mouse tissues using bisulfite sequencing. Our approach was powered to detect small RNA pseudouridylation based on robust detection of known {Psi} positions in tRNA fragments (tRFs), some of which show tissue-specific patterns of modification. In contrast with tRFs, we find that miRNA pseudouridylation is exceedingly rare, with a single modified miRNA (miR-3068-5p) identified in mouse tissues. Pseudouridylated miR-3068-5p diSerentially repressed predicted miRNA targets with less stable miRNA:mRNA pairing modes. This study fills a long-standing gap in transcriptome-wide {Psi} profiling and reveals a new potential function for {Psi} as a modulator of activity of small regulatory RNAs.

2
Pseudouridylation landscape across 42 S. cerevisiae cytosolic tRNA isoacceptors via Nanopore direct RNA sequencing

Barry, M. L.; Abu-Shumays, R. L.; Barnes, L. E.; Shaw, E. A.; Reinsch, J. L.; Vaaler, A. L.; Basham, Z. D.; Jain, M.; Koutmou, K. S.; Garcia, D. M.

2026-05-01 molecular biology 10.64898/2026.04.28.721490 medRxiv
Top 0.1%
19.1%
Show abstract

Pseudouridine is the most abundant RNA base modification due to its prevalence in tRNA and rRNA, where it serves as a key modulator of structure and function. Yet even in a widely used model organism, the budding yeast Saccharomyces cerevisiae, the positions of all pseudouridines in tRNA have not been completely annotated. Using Nanopore direct RNA sequencing (DRS), an established method for detecting RNA pseudouridylation positions, we sequenced cytosolic tRNA from eight pseudouridine synthase (PUS) knockout S. cerevisiae strains, including deletion strains of Pus1, Pus3, and Pus7. Analysis of these data verified thirty-four existing pseudouridine annotations and predicted eleven previously unannotated pseudouridine sites. Our analysis revealed DRS signal changes at several non-uridine sites with the loss of a PUS, including apparent changes in modification abundances at position 37 upon deletion of Pus3. LC-MS/MS and primer extension assays, however, indicated no change in the abundance of these modifications with the loss of Pus3. Our analysis underscores the need for caution in interpreting DRS-based signal changes, particularly in modification-dense regions. Combining existing modification annotations for the thirty-one isoacceptors in the Modomics database with our dataset that added annotations for the remaining eleven isoacceptors, we created a map of all detected pseudouridines, and the enzymes responsible for their catalysis, across the forty-two S. cerevisiae cytosolic tRNA isoacceptors.

3
The Illumina Stranded mRNA protocol is not strongly stranded for mRNA with low U content

Menshikova, O.; Nuez, I.; Courtier-Orgogozo, V.

2026-07-13 genomics 10.64898/2026.07.09.737636 medRxiv
Top 0.1%
10.8%
Show abstract

The Illumina TruSeq Stranded and Illumina Stranded mRNA protocols are commonly used for strand-specific bulk RNA-seq and they typically yield >99% antisense reads. We show that these protocols can generate sense-oriented reads for transcripts with extremely low U content (<3%). Indeed, such regions can bypass the dUTP-based blockade of cDNA second strand amplification. A small number of genes are affected by this issue (three in Drosophila melanogaster, including the glue gene Sgs3, and 46 in Mus musculus). To prevent overestimation of expression levels, we recommend excluding sense reads for all genes.

4
Systematic engineering and machine learning analysis of intrinsic terminators reveal crucial nucleotides directly upstream of the terminator hairpin.

Koster, C. C.; Terlouw, B.; Nieuwkoop, T.; Creutzburg, S. C. A.; Martin-Pascual, M.; Paredes Barrada, M.; Kopsiaftis, P.; Heilig, H. G. H. J.; van Laar, T.; van der Oost, J.; Claassens, N. J.

2026-07-07 molecular biology 10.64898/2026.07.06.736697 medRxiv
Top 0.1%
10.0%
Show abstract

Transcriptional termination efficiency is considered an important parameter for fine tuning bacterial gene expression. Still, the design principles that determine transcription termination efficiency remain poorly understood. In this study, we aimed to investigate the impact of the 3' untranslated region (3'UTR) on gene expression in Escherichia coli and other bacteria. First, 3'UTR variant sequences were generated, with randomized 30 bp sequences inserted between the STOP-codon and an intrinsic terminator, consisting of a GC-rich hairpin and a downstream poly(U)-tail. Using three reporter genes, it was found that different 3'UTR sequences resulted in an up to five-fold difference in protein production, independent of the upstream coding sequence. The highest protein production was achieved when an adenosine was present directly upstream of the terminator hairpin. This was consolidated by systematic substitution of key nucleotides of the terminator and assessing their effect on mRNA and protein levels. Subsequently, we developed a predictive random forest machine learning model trained on the termination efficiency of different natural and synthetic terminator sequences, revealing an important role for the nucleotides directly upstream of the terminator hairpin. Altogether, this study showed that an additional adenosine nucleotide upstream of the terminator hairpin leads to improved protein production while reducing terminator read-through.

5
Predicted Bacterial uRBSs Reveal Translational Coupling and Ribosome Mediated RBS Occlusion as Gene-Controlling Mechanisms

Dietz, T.; Hahnfeld, J. M.; Neumann, S.; Reinsch, Y. A.; Wenz, T.; Barth-Weber, S.; Blom, J.; Goesmann, A.; Evguenieva-Hackenberg, E.

2026-05-14 molecular biology 10.64898/2026.05.12.723736 medRxiv
Top 0.1%
9.7%
Show abstract

Upstream open reading frames (uORFs) in the 5' leader of bacterial mRNAs can modulate gene expression, yet genome-wide identification remains limited. We combined bioinformatic prediction of ribosome-binding sites (RBSs) - a Shine-Dalgarno sequence and a start codon - with experimental validation to uncover new uORFs in Sinorhizobium meliloti 2011. From totally 1106 predicted upstream RBSs (uRBSs), we first examined 15 candidates using eGFP reporters and integrating existing RNA-seq and Ribo-seq data. Translation was detected at 13 sites, with fluorescence intensity broadly correlating with predicted initiation rates. Two uRBSs correspond to gene start sites, thereby refining gene annotations. In nine cases, uRBS mutations affected downstream gene expression in reporter fusions. Among others, the data suggests that a Type I secretion system operon, the RNA chaperone gene hfq, and metabolic genes are regulated by uORFs. Four uORFs acted through translational coupling. We also identified uRBSs that were ribosome-occupied yet (nearly) silent in eGFP assays, and closely spaced to the downstream main RBS (mRBS). These uRBSs probably mediate ribosomal occlusion downregulating lacR and SM2011_RS36230. A re-screen of the prediction set revealed 335 close uRBS/mRBS pairs. Three of them were analyzed, supporting the proposed ribosomal occlusion mechanism for SM2011_RS03630 and SM2011_RS22110, while for glnK translational coupling to an uORF was suggested. These results indicate that uORFs are more widespread in bacteria than previously recognized and suggest that direct ribosomal occlusion of the mRBS is a novel mechanism for down-regulating protein synthesis.

6
Unexpected ribosome turnover during prolonged translation inhibition

Russell, P. J.; Clark, C. A.; Ashriem, M.; Kearse, M. G.

2026-05-07 molecular biology 10.64898/2026.05.06.723260 medRxiv
Top 0.1%
8.9%
Show abstract

Eukaryotes use several distinct quality control pathways to resolve aberrant ribosomes and mRNAs. For example, the no-go decay mRNA pathway is stimulated after ribosome collisions caused by stalled ribosomes translating damaged or truncated mRNAs. Separate decay pathways for non-functional 40S and 60S subunits containing rRNA mutations affecting decoding and peptidyl transferase activity, respectively, have also been elucidated. To our knowledge, whether eukaryotes have evolved a quality control pathway to sense and process globally stalled ribosomes is unclear; however, such a pathway would be advantageous to eukaryotes during exposure to natural elongation inhibitors such as ricin and diphtheria toxin. Here, we test how prolonged robust inhibition of elongation using a high dose of cycloheximide (CHX) affects ribosome turnover. Despite no decrease in cell viability and that mammalian ribosomes have been classically characterized of having a half-life of 3-5 days, a single 24 hr high dose of CHX resulted in drastically shortened half-lives (<24 hr) of 28S and 18S rRNA in A549 cells. A [~]2-fold reduction in nearly all ribosome species was observed by polysome analysis in HeLa and A549 cells after prolonged CHX treatment. Depletion of ribosomes was also evident when assessing ribosomal proteins from both the 40S and 60S subunits by Western blot. Literature supports that ribosomes can be degraded by autophagy and the ubiquitin (Ub)-proteasome system. Upon testing inhibitors of both pathways, only proteasome inhibitors (i.e., MG132 and bortezomib) rescued both rRNA and ribosomal protein levels. Proteasome inhibitors also rescued ribosome levels in polysome profiling experiments. Remarkably, rRNA levels were not rescued during CHX treatment when co-treated with the Ub activating enzyme E1 inhibitor, TAK243. Polysome analysis also showed that the high prolonged dose of CHX did not cause robust accumulation of collided ribosomes compared to control treatments. Proteasome-dependent turnover of rRNA was also observed with high doses of other elongation inhibitors, namely anisomycin, homoharringtonine, and lactimidomycin. The recognition capabilities of the pathway were further expanded as we observed that 80S ribosomes not trapped on the mRNA were also targeted for degradation by the proteasome. Together, our findings define the framework of a regulatory pathway in mammalian cells that degrades both ribosomal subunits in response to prolonged periods of robust elongation inhibition.

7
A pipeline for identifying small noncoding RNA (sRNA) candidates in bacteria

Elhedi, S.; NDiaye, K. D. S.; Perreault, J.

2026-07-10 molecular biology 10.64898/2026.07.02.735529 medRxiv
Top 0.1%
7.8%
Show abstract

Bacterial small non-coding RNAs (sRNAs) are central post-transcriptional regulators, yet their computational identification suffers from high false-positive rates due to transcriptional noise and the absence of canonical coding features. We developed a three-stage pipeline integrating sRNA prediction (sRNA-Detect), transcription start site mapping (TSSAR, dRNA-seq), and Rho-independent terminator detection (RNIE), applied across nine phylogenetically diverse bacterial species spanning six phyla. Sequential filtering achieved 1.4 to 33 fold precision improvements across nine species, reducing candidate sets by up to 99.6% while recovering known sRNAs at rates reflecting reference database depth (6% recall in S. aureus, 33-34% in E. coli and S. enterica) TSS and RIT constraints constitute universal, genome-size-independent biological filters that substantially enrich sRNA predictions across bacterial diversity. Precision variation across species reflects database incompleteness rather than pipeline failure, with unmatched predictions in poorly annotated organisms representing candidate novel sRNAs rather than false positives. RNA-seq coverage depth provides a reliable secondary indicator of biological relevance, though its interpretation requires accounting for sequencing depth variation across datasets.

8
Studying the regulons of OmrA and OmrB paralogous small RNAs reveals targets involved in central carbon metabolism and lipogenesis

Korepanov, A.;Jagodnik, J.;Quenette, F.;LAM, T.;HAMON, M.;Fromont, J.;Sismeiro, O.;Gherdol-Nouvion, V.;Maes, A.;Guillier, M.

2026-06-27 Molecular Biology 10.64898/2026.06.26.734639 medRxiv
Top 0.1%
7.7%
Show abstract

Small regulatory RNAs (sRNAs) are key players in bacterial adaptation to stress. They often occupy central positions in regulatory networks and control the expression of multiple targets. In a striking example of this, the enterobacterial OmrA and OmrB paralogous sRNAs are known to regulate about ten different targets, with extensive data suggesting the regulon is in fact much larger. Here we performed transcriptome and proteome analyses and identified more than fifteen new targets of Escherichia coli OmrA and OmrB. We validated several, including genes involved in central carbon metabolism and fatty acid synthesis, among which ppc, actP and fabA. Consistent with a role in carbon metabolism, overproducing OmrA or OmrB inhibited growth on glucose minimal medium. The analysis of suppressor mutants shows that this is due to a decreased carbon flux through the TCA cycle. Incorporating other datasets such as RIL-seq, we generated a multi-omics-based prediction of target candidates. Together, our results show that OmrA/B base-pair to various regions of their mRNA targets, and therefore likely act through diverse regulatory mechanisms. Hence, this work extends the OmrA and OmrB regulons, establishes an unsuspected connection with carbon usage, and shows the benefits of combining global analyses to investigate sRNA regulons.

9
OmrA sRNA Inhibits Translation of Phosphoenolpyruvate Carboxylase to Impair TCA-Cycle Flux

Stenum, T.; Le Huyen, K. B.; Kjellin, J.; Koskiniemi, S.; Wagner, E. G. H.; Holmqvist, E.

2026-06-26 microbiology 10.64898/2026.06.26.734723 medRxiv
Top 0.1%
7.4%
Show abstract

Small RNAs (sRNAs) rarely cause strong growth phenotypes upon overexpression, complicating efforts to link regulatory interactions to physiological outcomes. Here, we report that high levels of the Escherichia coli sRNA OmrA, but not its sibling OmrB, severely inhibit growth in glucose minimal medium. Genetic, biochemical, and physiological analyses indicate that OmrA-dependent toxicity results from reduced flux through the tricarboxylic acid (TCA) cycle. A UV-based suppressor screen identified mutations in the gene encoding Hfq, the RNA-chaperone that aids sRNA-mRNA interactions. Secondly, three independent mutations clustered in the ribosome-binding site of ppc, encoding phosphoenolpyruvate carboxylase, a key anaplerotic enzyme. OmrA directly inhibits Ppc translation via Hfq-dependent base-pairing in the ppc 5' UTR, including the mutated nucleotides obtained in the genetic screen. OmrA is significantly more effective than OmrB in ppc repression in vivo and in vitro, consistent with sequence divergence in their central regions. Supplementation with glutamate, glutamine, or downstream TCA cycle metabolites fully restores growth, linking reduced Ppc levels to metabolic limitation. These results identify ppc as a physiologically relevant OmrA target and suggest how RNA toxicity can uncover central metabolic nodes used by sRNAs to modulate bacterial physiology.

10
Translation of a four-codon ORF in the 3'UTR of MAP3K10 regulates its expression

Vasu, K.; Ghosh, S.; Sahoo, S.; Manna, D.; Som, S.; Akhtar, M. N.; Kar, D.; Eswarappa, S.

2026-04-29 molecular biology 10.64898/2026.04.27.721023 medRxiv
Top 0.1%
7.2%
Show abstract

Several eukaryotic mRNAs are polycistronic because of translatable open reading frames (ORFs) embedded in their untranslated regions (UTRs), in addition to the primary ORF (i.e., the coding sequence). While 5'UTR upstream ORFs (uORFs) are mechanistically and functionally well studied, 3'UTR downstream ORFs (dORFs) remain poorly understood. Here, we identify and characterize a short, translatable dORF in the 3'UTR of MAP3K10, which encodes a serine/threonine kinase involved in JNK signalling. A stringent computational screen predicted a conserved RNA G-quadruplex (rG4) within the 3'UTR of MAP3K10. Biophysical assays provided more evidence for rG4 formation, which drives translation of a conserved four-codon dORF. Disruption of the rG4 by point mutations or by an rG4-binding ligand reduced dORF translation. Reporter assays using constructs with strategically placed hairpin structures show that translation of the dORF is independent of both the 5 cap and the translation of the canonical ORF. Notably, deletion of the rG4-dORF module, either in exogenous constructs or in CRISPR-edited cells, led to reduced MAP3K10 expression. Together, these results provide evidence for the regulation of MAP3K10 expression by rG4-driven translation of a dORF. Thus, our study contributes to the growing body of evidence suggesting that dORF translation can regulate gene expression.

11
Dramatically reduced spliceosome, intronome, and splicing efficiency in Cyanidiococcus yangmingshanensis and Cyanidium caldarium

Slat, V. A.; Stark, M. R.; Rader, S.

2026-05-21 molecular biology 10.64898/2026.05.17.725761 medRxiv
Top 0.2%
6.2%
Show abstract

Eukaryotic pre-mRNA splicing is catalyzed by the spliceosome, whose ribonucleoprotein composition and the number of intron substrates it acts upon vary widely across eukaryotic lineages. The red alga Cyanidioschyzon merolae possesses a reduced spliceosome lacking the U1 snRNP, and an unusually small intron repertoire. We asked whether these traits are unique to C. merolae or shared across the related Cyanidiales and Cyanidioschyzonales lineages, as well as how they relate to splicing efficiency under light conditions relevant to photosynthetic growth. Genomic and transcriptomic analysis of C. merolae, Cyanidiococcus yangmingshanensis, and Cyanidium caldarium reveal that all three species harbour a reduced, but broadly conserved, set of splicing proteins. Strikingly, covariance model searches failed to detect U1 snRNA in either C. yangmingshanensis or C. caldarium, establishing U1 loss as a shared feature of all three lineages. We identified only 39 introns in C. merolae, 40 in C. yangmingshanensis, and 54 in C. caldarium. Splicing efficiencies were 42-50%, substantially lower than most organisms in which splicing has been measured, but low splicing is compensated by 2-4x higher expression of intron-containing genes than intron-lacking genes. Notably, light can enhance splicing efficiency in C. merolae and C. yangmingshanensis by up to 100%. Furthermore, the splice site and branch site consensus sequences are highly conserved and similar to those found in hemiascomycetous yeasts such as Saccharomyces cerevisiae. 85% of introns contain an in-frame stop codon with a strong bias towards the 5' end of the intron. These results indicate that dramatic streamlining of the spliceosome and intronome, together with inefficient splicing, predated the divergence of these lineages [~]320 million years ago, and is therefore a defining molecular trait of these extremophilic red algae.

12
RT-nested and interfering-Primer PCR reveal prevalent isoform-specific A-to-I RNA editing in neuronal genes

Wang, Z.; Ni, Y.; Cai, W.; Li, H.; Duan, Y.

2026-05-17 molecular biology 10.64898/2026.05.15.725286 medRxiv
Top 0.2%
5.7%
Show abstract

BackgroundMetazoan adenosine-to-inosine (A-to-I) mRNA editing temporospatially diversifies the neuronal transcriptome and proteome. The limited read length from next-generation sequencing (NGS) constrains the quantification of the potentially differential editing levels across different splicing isoforms, restricting our understanding of the extent to which RNA editing contributes to molecular diversity and its interplay with splicing. MethodsWe employed reverse transcription nested PCR (RT-nPCR) and developed a novel interfering-Primer PCR (iPrimer PCR) technique to distinguish different transcripts of any gene. We selected multiple essential genes exhibiting RNA editing in coding sequences (CDSs) or untranslated regions (UTRs) for isoform-specific amplification and Sanger sequencing. ResultsNine different Adar isoforms together with pre-mRNA had distinct editing levels at the S>G auto-recoding site, which was predicted to have isoform-specific effects on catalytic activities. Although pre-mRNA editing might exert isoform-dependent promotion/suppression of splicing, closely located editing sites, such as those in neuronal genes qvr and stj, still exhibited high correlation in editing levels due to co-editing. iPrimer strategy further discovered differential recoding levels between the long/short 3UTR isoforms of gene jef. ConclusionsWe provide the first comprehensive solution for isoform-specific PCR amplification of any gene, enabling quantification of RNA editing level of different isoforms. Our results offer insights into how RNA editing interplays with splicing, and highlight its complicated role in expanding molecular diversity. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=79 SRC="FIGDIR/small/725286v1_ufig1.gif" ALT="Figure 1"> View larger version (17K): org.highwire.dtl.DTLVardef@1ebc82org.highwire.dtl.DTLVardef@1ea365dorg.highwire.dtl.DTLVardef@1971aceorg.highwire.dtl.DTLVardef@160d053_HPS_FORMAT_FIGEXP M_FIG C_FIG We developed isoform-specific PCR followed by Sanger sequencing, and achieved the quantification of differential RNA editing levels in different transcripts of a gene.

13
Identification of implications of m6A regulators and autophagy-associated genes for prognosis in ovarian cancer

Chen, Y.; Yu, X.; Chu, W.; Shang, S.; He, N.; guo, l.

2026-06-29 obstetrics and gynecology 10.64898/2026.06.25.26356535 medRxiv
Top 0.2%
5.6%
Show abstract

The most prevalent RNA alteration in the mammalian genome is N-6-methylenediosine (m6A). There is mounting evidence linking dysregulation of m6A regulatory factors and alterations in m6A levels to the development, course, or prognosis of ovarian cancer. Genes having prognostic value were screened using the univariate, multifactorial, and Least Absolute Shrinkage Selection Operator (LASSO) Cox regression analyses. Important genes' m6A expression in clinical material was verified by real-time fluorescent quantitative polymerase chain reaction (RT-qPCR). In present study, all 23 regulators were significantly differentially expressed in ovarian cancer tissues. LASSO regression analysis screened for 10 key genes associ-ated with both autophagy and m6A. A risk score was constructed and nomogram was developed to forecast the prognosis of ovarian cancer patients. Additionally, individuals with ovarian cancer were classified as high-risk or low-risk; and the low-risk group might be more likely to benefit from im-munotherapy. RT-qPCR was used for the bioinformatics study of human ovarian cancer and normal tissues. Lastly, PLK2 and LEPR were confirmed to be associated with tumorigenesis in scRNA-seq. The risk score established by m6A and autophagy can be used to predict prognosis and susceptibility to anticancer drugs in patients with ovarian cancer.

14
Genetic detection of RNA-protein interactions using a bacterial three-hybrid assay

Gravel, C. M.; Berry, K. E.

2026-07-09 molecular biology 10.64898/2026.06.26.734845 medRxiv
Top 0.2%
5.5%
Show abstract

The bacterial three-hybrid (B3H) assay is a powerful genetic tool for detecting interactions between RNA and RNA-binding proteins (RBPs) and assessing the consequences of RBP mutations. This transcription-based system connects the strength of an RNA-protein interaction to the expression of a lacZ reporter gene in Escherichia coli cells. This in vivo approach allows researchers to dissect RNA-protein interactions within a cellular environment, bypassing the need for biochemical purification of RNAs or proteins. This chapter details a three-day protocol for generating quantitative B3H data. Since a significant challenge in B3H assays is RNA misfolding, we describe a recently optimized set of B3H constructs that mitigates this issue by isolating bait RNAs as discrete folding units.

15
Characterisation of the RNA-Binding Properties of the MRSA β-lactam resistance enzyme PBP2a

Christopoulou, N.; Dương, N. H.; Arede-Rei, P.; Torrens, G.; Blandenet, M.; Cava, F.; Granneman, S.

2026-07-07 biochemistry 10.64898/2026.07.05.736576 medRxiv
Top 0.2%
5.4%
Show abstract

Analysis of RNA-binding proteome data from different bacterial species revealed many cell wall metabolic enzymes cross-linking to RNA in vivo, hinting that these proteins directly bind RNA. Surprisingly, penicillin-binding proteins (PBPs) were also abundantly identified as putative RNA-binding proteins. The cell surface localisation properties of many of these proteins therefore beg the question at what stage of their cellular life cycle these proteins interact with RNA and what the functional significance is. Here, we characterised the RNA-binding activity of PBP2a, the alternative transpeptidase that confers {beta}-lactam resistance in MRSA. Using in vivo RNA-binding assays, we show that PBP2a interacts with hundreds of transcripts without apparent sequence specificity. Computational analyses identified a possible RNA-binding cleft in PBP2a proximal to its active site. Mutation of only two predicted positively charged residues located in this cleft substantially reduced cross-linking in vivo, implying that RNA recognition is largely dictated by RNA backbone interactions. While PBP2a does not regulate RNA steady-state levels, RNA-binding appears important for proper protein function: an RNA-binding deficient mutant exhibits reduced oxacillin resistance. These findings establish PBP2a as an RNA-binding protein in vivo and provide a framework to investigate how this non-canonical interaction may relate to cell wall biogenesis and {beta}-lactam resistance.

16
Estimation of splicing metrics for NMD-sensitive transcripts

Zavileyskiy, L.; Vlasenok, M.; Kuznetsova, A.; Skvortsov, D. A.; Pervouchine, D. D.

2026-07-04 bioinformatics 10.64898/2026.06.30.735642 medRxiv
Top 0.2%
5.3%
Show abstract

Alternative splicing is commonly quantified using the Percent-Spliced-In (PSI) metric, which measures the relative abundances of alternatively spliced isoforms. However, some transcript isoforms are targeted by the nonsense-mediated decay (NMD) pathway, introducing a strong bias that leads to underestimation of their true splicing rates. To correct for this bias, we developed an analytical framework and a set of statistical models employing a linear fractional transformation depending on a single parameter capturing the degradation rate of NMD-sensitive transcripts relative to normal mRNA decay. Using Gaussian mixture models, we demonstrated a clear separation of splicing events into two classes, responders and non-responders, with the former exhibiting strong upregulation upon NMD inhibition and the latter showing little or no response. Moreover, non-responders displayed higher coding potential and stronger translation signals both upstream and downstream of the stop codon, which are characteristic of NMD escape through translational readthrough. We further showed that incorporation of event-specific relative decay rates improves the interpretation of differential splicing patterns for NMD-sensitive transcripts. In sum, our results provide a solid framework for unbiased estimation of splicing metrics in NMD-sensitive transcripts from short-read RNA-seq data, without requiring NMD inhibition experiments.

17
Characterizing the Small Non-Coding RNA Pathways in the Invasive Zebra Mussel (Dreissena polymorpha)

Hernandez Elizarraga, V. H.; O'Brien, L. G.; Ballantyne, S.; Gohl, D. M.

2026-07-11 genomics 10.64898/2026.07.10.737777 medRxiv
Top 0.2%
4.4%
Show abstract

The zebra mussel (Dreissena polymorpha) is an invasive species that causes extensive economic and ecological damage. Here, we identify and characterize the key components of the small RNA (sRNA) and RNA interference (RNAi) pathways in zebra mussels. Like other mollusks, zebra mussels have extensive microRNA (miRNA) and Piwi-interacting RNA (piRNA) machinery but lack or have modified canonical factors needed to produce small interfering RNA (siRNA). Specifically, the zebra mussel Dicer sequence displays substitutions in the conserved DEAD box motif that is required for substrate processivity, and this organism also lacks some attendant accessory factors such as R2D2. We sequenced the small RNA found in both isolated somatic tissue (adductor muscle) and whole animals (including germline), and identified both conserved and novel miRNA and diverse piRNA sequences, but few endogenous siRNAs. To determine whether their remaining sRNA machinery could still be co-opted to initiate gene silencing, we injected dsRNA targeting several genes into zebra mussel adductor muscle. The injected rpn8-targeting dsRNA reduced rpn8 mRNA levels and was processed into sRNA that resemble endogenous miRNAs and piRNAs. The levels of both sRNA types correlated with mRNA knockdown, suggesting that they may act together to initiate RNAi as seen elsewhere. dsRNA targeting other genes produced variable results suggesting that particular criteria may be needed to trigger an RNAi response in this assay. Our results characterize endogenous sRNA pathways in zebra mussels, establish that dsRNA can induce RNAi, and lay the groundwork for further optimizations to establish RNAi-based genetic manipulation tools for this damaging invasive species.

18
misoTar: A novel approach for predicting miRNA and isomiR targets

Ripan, R. C.; Li, x.; Hu, H.

2026-05-12 bioinformatics 10.64898/2026.05.08.723919 medRxiv
Top 0.2%
4.3%
Show abstract

Understanding the interactions between microRNAs/isomiRs and mRNAs has long been a major challenge in RNA biology. Although numerous computational approaches have been developed to predict these interactions, most fail to account for isomiR mediated targeting. To address this limitation, we developed misoTar, a deep learning framework trained on more than 6.662 million positive and negative interaction pairs derived from 67 publicly available human samples across six independent studies. In five-fold cross-validation, misoTar achieved an average precision of 0.930 and a recall of 0.898. Evaluation on independent test datasets demonstrated consistently superior or comparable performance relative to existing tools, including TargetScan, Mimosa, DMISO, and TEC-miTarget. In addition, single-nucleotide mutation analyses of true positive interactions revealed the critical functional contributions of non-seed regions in microRNA/isomiR targeting. Overall, misoTar provides a robust and accurate framework for predicting microRNA/isomiR interactions while offering new biological insights into microRNA targeting mechanisms. The misoTar tool is publicly available at https://figshare.com/projects/misoTar/262723.

19
Identification of Potential Regulatory Non-Coding RNAs in Lotus Japonicus Symbiosis

Budnick, A.; Utley, D.; Blahovska, Z.; Radutoiu, S.; Sederoff, H.

2026-05-21 plant biology 10.64898/2026.05.19.726297 medRxiv
Top 0.3%
4.2%
Show abstract

O_LISymbiosis between legumes and rhizobia is beneficial on nutrient-poor soils, as it enables the fixation of atmospheric N2. To establish this symbiosis, gene expression in both the host plant and the symbiont has to be regulated. To understand the underlying RNA-mediated regulation of host gene expression, we designed experiments to identify competing endogenous networks involving circular RNA, microRNA, and linear transcripts during symbiosis, using wt and symbiosis-deficient Lotus japonicus mutants with the rhizobium Mesorhizobium loti (M. loti). C_LIO_LICircRNA, miRNA, and linear transcripts were identified from Lotus japonicus wildtype and CCamK mutant (ccamk-13; snf-1) seedlings without inoculation or with M. loti inoculation using deep short-read sequencing with rRNA-depletion and random primers. C_LIO_LIDifferentially expressed miRNAs showed negative correlations to predicted target genes and may regulate symbiotic processes. The symbiosis essential iron-sensor LjnsRING/BRUTUS expresses a circRNA which was upregulated in symbiotic treatments. This circRNA may act as a target mimic and contribute to nodule longevity. CircRNAs are predicted to act predominantly as trans-regulatory molecules with similar frequencies in Arabidopsis thaliania, Oryza sativa, and Lotus japonicus. C_LIO_LIWe identified novel miRNAs, long noncoding RNAs, and circRNAs, and nominated several as potential new regulatory non-coding RNAs that may act as target mimics to stabilize genes and support symbiosis. C_LI SummarySymbiosis between Lotus japonicus and Mesorhizobium loti involves treatment-specific regulation of competing endogenous RNA networks involving circular RNA, miRNA, and linear transcripts.

20
A Metabolic Enzyme, Pyruvate Carboxylase, Functions as a Sequence-Selective Small RNA Sensor for Antiviral Immunity

Kariyawasam, U.; Goswami, S.; Hao, M.; Wiscovitch-Russo, R.; Chen, Q.; Yang, J.; Qiu, J.; Marquez, M.; Sui, H.; Chang, W.; Imamichi, T.

2026-07-03 immunology 10.64898/2026.06.29.735367 medRxiv
Top 0.3%
4.1%
Show abstract

Interleukin 27 (IL-27) is an anti-HIV cytokine that induces 14 novel microRNAs (miRNAs) in T cells. We previously reported that transfection of two of these miRNA mimics, miRTC10 and miRTC14, differentially induced interferons (IFN)A2, A8, A13, and L1 expression in human primary macrophages. However, the mechanism underlying this activation remains unclear. Here, we show that miRTC14 does not directly target IFN-regulatory genes but instead engages cytosolic RNA-sensing proteins. Using miRNA pull-down coupled with mass spectrometry and immunoblotting, we identified a metabolic enzyme, pyruvate carboxylase (PC) and laboratory of genetics and physiology 2 (LGP2/DHX58) as direct binding partners of miRTC14. Functional analyses revealed that miRTC14 induces IFN expression by more than100-fold (p < 0.001), whereas PC and LGP2 depletion markedly attenuated this response (50-100 fold reduction, p < 0.01). Reconstitution of PC and LGP2 in deficient HEK293 cells restored miRTC14-driven IFN induction. We found that miRTC14-induced IFN activation depends on sequence features at the duplex termini and is unlikely to arise from canonical miRNA-mediated gene silencing. These findings establish PC as a novel miRNA-binding protein and define a previously unrecognized RNA-sensing mechanism by which miRTC14 drives IFN production, linking metabolic enzymes to RNA sequence-dependent innate immunity.